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The concepts, which have traditionally been useful in understanding the effects of the electron- 
phonon interaction in optical spectroscopy, are based on insights obtained within the infinite elec- 
tronic band approximation and no longer apply in finite band metals. Impurity and phonon con- 
tributions to electron scattering are not additive and the apparent strength of the coupling to the 
phonon degrees of freedom is substantially reduced with increased elastic scattering. The optical 
mass renormalization changes sign with increasing frequency and the optical scattering rate never 
reaches its high frequency quasiparticle value which itself is also reduced below its infinite band 
value. 
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I. INTRODUCTION 

Many of the physical insights that have guided the interpretation of data on effects of the electron-phonon interaction 
in metals are based on an infinite band model with a constant featureless electronic density of states (EDOS)^. fn the 
eighties there appeared several studies, mainly motivated by the physics of the A15 compounds, which took account of 
energy dependence in the EDOS around the chemical potentia l^i^i"^ . Band structure calculations for the A15 showed 
peaks in the EDOS with variations on the energy scale of 50 meV. A variety of experiments also showed sensitivity 
of properties to disorder. For example, disordered MoaGe has a higher value of superconducting critical temperature 
than its crystalline counterpart. This is naturally explained if the chemical potential in ordered metallic MosGe falls 
in a valley of the EDOS. Radiation damage then fills this valley and leads to an increase in EDOS at the Fermi energy 
and a higher value of T^. 

More recently several authors have considered a different but related effect, namely a finite band^^SiL^. One 
experimental realization of this situation is the fuUeride compounds MaCgo {M - an alkali metal), where band 
structure calculations^ show narrow band with the width W of the order of 1 eV, while the phonon spectrum extends 
up to about 200 meV in some cases. Physical consequences brought about by the finite band can be studied in the 
frame of a simplified (particle-hole symmetric) model with a constant Nq with a cut off applied at zLW/2 where W 
the band width related to A''o by A'o = 1/W. A somewhat surprising result of such studies is that, even for rather 
wide bands {W of order a few eV) certain aspects associated with the effect of the electron-phonon interaction are 
profoundly modified as compared to the corresponding infinite band behavior. For example, in an infinite band with 
constant electronic density of states, the electron-phonon interaction leaves iVo(w) unaltered and no phonon structure 
appears in the dressed normal state EDOS. To see phonon structure it is necessary to go to the superconducting state 
which develops a gap and consequently a non constant EDOS. However if a cut off is applied to the constant Nq{uj), 
then phonon structure appears in the dressed quasiparticle density of states as it does in the superconducting state 
and also in any case when EDOS is non constant around the Fermi energy. The phonon structure which appears 
in the dressed EDOS is surprisingly significant in magnitude even for modest value of the electron-phonon mass 
renormalization parameter A. Mathematically the self energy must be solved for self consistently when a finite band 
cut off is introduced. This contrasts with the infinite band case where the bare Green's function can be employed in 
the self energy expression--. Self consistency leads to a smearing of the band edge region as well as a widening of the 
band. As the total number of states in the electronic density of state must remain constant, this transfer of spectral 
weight to higher energies beyond the bare EDOS cut off, implies that it must correspondingly be reduced at smaller 
energies and the details of this reduction depend significantly on the phonon energy scale and coupling strength to 
the various phonon modes. The effect of widening of a finite electronic band due to the electron-phonon interaction 
has been observed and discussed previously by Liechtenstein et aiii in the context of fuUeride compounds. 

There are many qualitative changes in electron-phonon renormalization effects which have their origin in finite 
bands. For example the real part of the electronic self energy Si(u;) for w > is everywhere negative in an infinite 
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band and decays to zero beyond a few times the maximum phonon energy, which we denote cj/j. Thus the electronic 
effective mass is always increased by the electron-phonon interaction and returns to its bare mass value from above 
at a few times ujd- By contrast for a finite band, as described in Ref. change sign as uj increases and the 

renormalized mass at high uj can actually be smaller than the bare band mass. This is an example of qualitative 
change brought about in the electronic self energy by finite band effects. Others are described in the recent paper of 
Cappelluti and Pictronero^, who also considered the effect of impurities. 

In this paper we consider optical properties with particular emphasis on the combined effect of temperature and 
impurity scattering in a finite band electron-phonon system. In Section II we provide a brief summary of the formalism 
needed to compute the electron self energy S vs w for a system of electrons coupled both to phonons and to impurities. 
We also present analytic formulas which apply in the non selfconsistent approximation. They will prove useful for 
interpretation of the numerical results. The optical conductivity without vertex corrections follows from the Kubo 
formula for the current-current correlation function. The optical self energy, or the memory function, is introduced 
and related to the complex optical conductivity cr(ti;). We summarize some known approximate but analytic formulas 
for the optical scattering rate and effective mass renormalization which have been found useful in past studies related to 
infinite (very wide) electronic bands. These formulas appropriately modified in the context of finite bands are applied 
to obtain a description of the non selfconsistent approximation, which provide insight into the various features found in 
numerical solution of the full equations. Two models for the electron-phonon spectral density a^i^(w) are introduced. 
For definiteness both are based on the specific phonon spectrum of K^Cqq. One consists of three delta functions 
suitably chosen to mimic the real spectrum while the other one uses truncated Lorentzians instead of delta functions 
to help understand the modification brought about when the extended nature of real spectra is accounted for. In 
Section III we describe results for the case of a rather wide band and the three delta function model for a^F(w) with 
a modest value of A = 0.71. We start with a discussion of the dressed electronic density of states with emphasis on 
temperature and impurity effects. Then the memory function is analyzed and compared with the quasiparticle self 
energy, the differences arising from finite band effects are emphasized. In Section IV we present the results for an 
extended electron-phonon spectrum, increasing the spectral A and decreasing the width of the band. Section V is our 
conclusions. 

II. FORMALISM 

A. The electronic self energy and the renormalized EDOS 

The central quantity of our problem is the electronic self energy ^(z) = Si(z) + 1^2(2) ■ It is calculated from the 
Migdal equations formulated in the mixed real-imaginary axis representatioii^*i^: 
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where LUm = 7rT(2m — 1), m Cz Z are the fermionic Matsubara frequencies, and f{Lo) and n{uj) are the Fermi and 
Bose distribution functions respectively. The electron-phonon interaction is specified in terms of the electron-phonon 
spectral function a^F(uj) (the Eliashberg function). The parameter F, which has the meaning of an impurity scattering 
rate, specifies the strength of the interaction with impurities. The variable z in eqs. can assume arbitrary 

complex values. Description of spectroscopic experiments requires knowledge of the retarded electronic self energy 
at real frequencies, which corresponds to solutions with z = lo -\- iO^. A fast and stable numerical procedure for this 
purpose was proposed by Marsigiio et afl^. It starts with computing the solutions for Tt(z) on the imaginary axis, 
at z = iujmi where eq. ^ is simpler. Then, the function rjliujm) is used to set up an iterative procedure to find 
T,{iLj + iO"*") just above the real axis. 

The quantity Nq{^) appearing in Eq. © is the bare EDOS. In this paper we use for it the following simple model: 



NoiO = N,e {W/2 - 1^1) , 



(4) 
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where W is the bare band width and 8(x) is the step function. The constant A''o is fixed by normahzation: Nq = 1/W. 
In this paper we retain particle-hole symmetry for simplicity, with the chemical potential at the center of the band, 
/i = 0. In the clean case it was showni that the main characteristic features appearing in the electronic self energy 
and the memory function due to the finite band width, do not depend significantly on details of the bare electronic 
band. 

The renormalized density of electronic states, or density of states for quasiparticles, is defined by 

/ + 00 
dCNo{OA{^,iu). (5) 
-OO 

Here 



A{^,Lo) = -lmGret{Lu;)/7r (6) 
is the electronic spectral density and the retarded Green's function Gret{£,T^) is defined by the relation 



[Gret{.L0j)V' ^{[G^{i,z)]-' -nz)] (7) 

with Go{^,z) = l/{z — ^) being the free electron Green's function. The renormalized quasiparticle can be expressed 
in terms of the function rj — rji + ir]2 of Eq. ||2J) as follows N{u!)/No{0) — ~ri2{uj)/TT. 

The renormalized density of states N{uj) is a very important quantity. It features various signatures of the interaction 
of electrons with phonons and impurities. Note that in the infinite electronic band approximation with flat bare EDOS 
the renormalized EDOS N(uj) remains constant and does not carry any physical information. In the present context 
of a finite band, N{llj) for electron-phonon system has been studied recently by Dogan and Marsiglic^ at T = 
and by Knigavko and Carbottai at finite temperatures. Below we emphasize the analysis of the combined effect of 
both phonons and impurities. We find that knowledge of the features of the renormalized EDOS helps to understand 
better the behavior of the other spectroscopic quantities such as the memory functions, which are related to the 
optical response. The renormalized EDOS N{uj) itself is a measurable quantity and can be directly probed by 
tunneling spectroscopy or angle-integrated photoemission spectroscop}»i2iiiiii. The accuracy of the latter technique 
has increased dramatically in recent years and properties of both new and traditional materials have been scrutinized. 
It has been argued in Ref. 8 that normal state boson structure should be detectable in such experiments for metals 
with electronic band width of order a few eV. 

Let us return to the electronic self energy. For the purpose of the following discussion we present the general Eq. 
([T|l for the values of the argument just above the real axis, namely z ~ uj + iO+ and at temperature T = (note that 
henceforth we will use real axis variable, such as w, as shorthand for uj + iO^). Separating the real and imaginary 
parts of E we obtain the following expressions: 

El w = TP \-4 + 2luP dJ \ [ / dn— , (8 

^' J.^LJ-u;' No{0) Jo No{0)Jo lj^ ~ {uj' + Q)^ ' 



(9) 



where the symbol P in Eq. ||SJ) means that in the divergent integrals over lu' the Cauchy principal value has to be 
taken. Note that the self consistent nature of this equations is now masked. On the right hand side of Eqs. lO and 
© the self energy enters only via the renormalized EDOS N{uj). 
The quasiparticle mass renormalization is defined by the relation: 



X(fff) = - hm dEi(tj)/da;. (10) 



From Eq. we find that it is given by 



J^^ CO No{0) Jo ^ 'Jo (cu + nrNoiOy ^ ' 

where N'{uj) = dN{uj)/doj is the derivative of the renormalized EDOS. For an infinite band with a constant N{lu) we 
recover the known result. The second term on the right hand side of Eq. (|ll|l reduces to A = 2 dfla^ F{fl) /fl, the 
usual expression for the mass renormalization due to electron-phonon interaction, while the first term, which represents 
the effect of the elastic scattering from impurities, vanishes. In the case of an energy dependent EDOS impurities 
produce a finite contribution to the quasiparticle mass renormalization. Our subsequent numerical analysis shows that 
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for a finite band N'{u!) is a decreasing function of w for w > in large intervals, which become especially substantial 
a u!d ^ W [remember that N{—u!) = N{uj) because we consider the half filling case]. This makes the impurity 

contribution to X^qJ^^ negative and opposite in sign to the phonon contribution. Therefore, in a finite electronic band 
the increased elastic scattering results in the apparent decrease of the magnitude of the electron-phonon interaction, 
as specified by Xq'jf^^ ■ 

In this paper we solve the self consistent equations for the self energy numerically. To better understand the trends 
observed in our numerical results it is helpful to have analytic, though maybe not exact, expressions for the self energy. 
We found that a useful approximation is to replace the renormalized EDOS N{oj) in Eq. ^ and Hll|) with the 
bare EDOS No{uj), given in the model we consider by a constant equal to A^o(O) with cutoff at the bare band edge 
W/2 (see Eq. (0}). This approximation amounts to disregarding the self consistency, while keeping track of the finite 
width of the band. It is expected, and we confirm this in the following Section, that this approximation is good at 
small frequencies uj as long as the characteristic phonon frequency is much smaller than the band width W. Moreover, 
this non selfconsistent approximation allows us to obtain a simple estimate for the characteristic frequency ujqp of a 
finite electronic band, when the real part of the self energy changes sign. Deficiencies of the non selfconsistent 

approximation are discussed later. The non selfconsistent results have the form: 



■V (ns) 

qp 



= Tin 



= -2 



Lu + W/2 



Lu - W/2 



W/2 + Q 



Te{W/2- \uj\) 



Jo 

/•LU 

-/ dna^F{n)e{w/2-\uj-n\ 

Jo 



uj + n UJ- w/2 - n 



w/2 



1 



Q l + n/(W/2)' 



(12) 
(13) 
(14) 



These equations, while not exact, reduce to the well known infinite band results in the limit W oo and show 
the modifications brought about by a finite band. Note in particular that when no impurities are present the mass 
renormalization is always reduced over its infinite band value and that the quasiparticlc scattering rate t~j^(lj) = 
— 2I]2(i-l' + iO+) drops to zero for uj > W/2 + lud instead of remaining constant as in the infinite band case (we remind 
that luj:) denotes the maximum phonon frequency in a^F{fl)). 



B. Optical response 



For characterization of the optical response the quantity of interest is the memory function, which is the optical 
counterpart of the self energy. The memory function M{u!) = Mi{ui) + iM2{uj) appears explicitly in the following 
expression for the complex optical conductivity (t(w) = cri{uj) + icr2(a;): 

o Q 1 

a{u;) = — . , (15) 

TT M (Uj) — tUJ 

which is also called in the literature the extended Drude formulaiS^iliiSiiSiSfi. In this equation S is the optical sum 
defined by the integral 

ai(Lu)duj. (16) 

For the infinite and flat electronic band it was shown by AllenSi that in the limit u ^ the memory function and the 
self energy are closely related, namely: Mi(0) — —2^2(0) and M2(0) = Si(0). This is one of the reasons why the real 
part of the memory function can be identified as the optical scattering rate, Mi{oj) = t~^{u)). On the other hand, the 
imaginary part of the memory function is usually related to the frequency dependent optical mass renormalization, 
M2{ijj) = —uj\op{ijj) ■ The memory function can easily be found if the conductivity is known. Indeed from Eq. 115|l it 
follows that 



, _ 25 o-i(aj) 

^ TT (Tl((-0) + 0-2(0;) 



\ ( \ 25 cr2(w) 

-wAop(cj) = UJ — J—- (18) 



and these are the relations that we used in our numerical work presented below. 
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The optical conductivity was obtained using linear response theory, neglecting vertex corrections. The details were 
described previously elsewhere^, and here we just write down the expression for the real part of the conductivity: 



f {u') - f{L0' + uo) 



(19) 



^ is the averaged over the Brillouin 
In the following discussion of the 



and point out that the corresponding imaginary part was computed as the Hilbert transform of the real part, based 
on the the Kramers-Kronig relations. In Eq. H19|) N^{£) = iVo(^)w| where v 
zone the square of the group velocity defined for a general dispersion in Ref. |2: 
optical response and in the numerical calculations of this paper we assume that the system is isotropic and use for 
the expression u| — (-^ + derived from the quadratic dispersion of free electrons with lower band edge at 
^ = —W/2 {D is the number of spatial dimensions, rn the free electron mass). It is useful to introduce the optical 
effective mass renormalization as 



0), 



(20) 



which is the quantity to be compared with the quasiparticle effective mass renormalization Xqp of Eq. (|10|l . We 
will see that complete numerical results indicate that these two quantities are not the same in a finite band. Another 
important quantity that we will discuss below is the optical scattering rate at the Fermi level, T~p{ijj = 0). 

As we have done for the self energy it is helpful in understanding the complete numerical results, that will be 
presented in the following sections, to have simple although approximate analytic expressions for the optical quantities 
with which to compare. It is not feasible to obtain a simple accurate expression for the optical scattering rate in the 
general case, and even our model of bare EDOS with sharp cutoffs (see Eq. Q) do not provide enough simplifications 
for this purpose. We decided therefore to make use of the existing expressions, valid for the infinite band. Historically, 
based on second order perturbation theory for the electron-phonon system Allen^'^ was the first to provide such an 
equation valid at zero temperature in the infinite flat band case. A generalization to finite temperature was made 
by Shulga et aU^ using a very different method which starts with the Kubo formula and makes approximations to 
get the same result as Allen when T ^ limit is taken. On the other hand Mitrovic and FiorucciS^ and Mitrovic 
and Perkowita2& have generalized Allen's original work to include the possibility of an energy dependent electronic 
density of states. They considered only zero temperature. Very recently Sharapov and Carbotta^i have provided a 
finite temperature extension based on the Kubo formula. Such formulas have also been used recently in analysis of 
datck^Si^ and in comparison with more complete approaches^ii^S. The formulas for the optical effective mass and 
scattering rate, which are suitable for our forthcoming discussion, are those for T = 0. They are given in Refs.l25l and 
and here we reproduce them for the reader's convenience: 
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where the symbol P means, as usual, that the lu' integrals are calculated as the principal Cauchy values. In the above 
equations N{llj) is the renormalized EDOS that fully incorporates the self consistent electronic self energy. To grasp 
the finite band effects in the memory function we intend to replace N{u!) with the bare EDOS Nq{llj) from Eq. Q 
similarly to our approach to the derivation of Eqs. (|12|I - H13|) for the non self consistent self energy. Note however that 
in order to arrive at Eqs. (|21|l and H22|l it is necessary to assume^'"' that iV„(^) is constant and extends to infinity, 
i. e. no cut off is applied to it. This means that after the proposed replacement effectively only finite band effects 
originating in the self energy are included in the optical quantities and the resulting approximate formulae are not 
expected to be quantitatively correct for all frequencies lo. Nevertheless we found that such analytical expressions are 
very useful at a; < W/2. They read: 



Aop(w) = — 



In 



W/2-UJ 



UJ 



W/2 + UJ 



w/2 



UJ 



In 



\WI2) 



UJ 



dna^F{n) 



In 



LU -Q- w/2 uj + n 



LU + Q + W/2 u-n 



Lu^~{n + w/2f 



{n + w/2f 



w/2 



In 



(SI + W/2f 



_ (r! + W/2f 



0(w) 



27rr 



e(VK/2 -uj) + ^9(0' - w/2) 



— / dn a^F{n)Q{Lu - n) [{oj - n)Q{w/2 - (w - n)) + T4^/2e(t^ -n- w/2)] 

^ Jo 



(23) 



(24) 



6 



In particular, Eq. H23|) is used below to obtain a reasonable estimate for the characteristic frequency ujop at which the 
imaginary part of the memory function changes sign. 



III. RESULTS DUE TO A BAND CUTOFF 



Motivated by the elect ron-phonon interaction in the fuUeride compound KaCgo, we use a three frequency model 
for the electron-phonon spectral function: 



^ 2 

1=1 



LO - 



(25) 



with X]^=i ^1 = 1- The interaction strength a is defined as the area under the c?F{iS) curve. The mass enhancement 
parameter A is given by eq. (0) with z = 0. We set A = 0.71 with li = 0.3, 12 = 0.2, Z3 = 0.5 and loi : 102 : lo^ = 
0.04 : 0.09 : 0.19 eV^. This model has a = 43.8 meV and win — 102.5 meV, where win is the logarithmic frequenc}«i, 
a convenient parameter to quantify the phonon energy scale. For the forthcoming discussion we set W = 2.5 eV, 
which leads to a small value for the adiabatic parameter win/(W^/2) = 0.082. For the model of Eq. 1)251 the non self 
consistent approximation to the self energy given by Eqs. 1)1211 - 1)1411 becomes 
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(26) 
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For the parameters chosen the non selfconsistent mass renormalization is 0.63 in the clean case to be compared with 
A = 0.71. 

We begin by reviewing effects of the electron-phonon interacti on d ue to a finite bandwidth which are seen in the 
EDOS. Some of the features have been studied previously in Refs. l6l7l H. Here we want to emphasize impurity effects 
and give the comparison between non selfconsistent and fully selfconsistent results. In Fig.n(top frame) we show the 
frequency dependence of the renormalized quasiparticle density of states N(uj) based on the three frequency model 
of Eq. ()25)l with A = 0.71 and a half band width W/2 = 1.25 eV. These parameters are by no means extreme yet 
the deviations from the infinite band case [N{u!) = 1 for all uj] are substantial. First, note that a three step phonon 
structure is clearly seen at small lu in the lower temperature curves. The top set of four curves (solid) are for F = 0, no 
residual scattering, and the four lower curves (dashed) are for F — 22.2 meV or a residual scattering rate of 140 meV 
for the chosen value of W. The temperatures are T = 14.5, 72.5, 145 and 435 K. For the 435 K curve the thermal 
smearing is large but not for the others. As the impurity scattering rate is increased the band width increases but 
the phonon structures do not smear appreciably. Instead their relative amplitude is slightly attenuated. Identifying 
the three low frequency plateaus in N(uj) we plot, in the bottom frame of Fig.^ their heights (solid stars, squares, 
triangles, refer to the left vertical axis) as a function of F and compare with the value of EDOS at a; = (solid 
diamonds, refer to the left vertical axis). All are reduced in magnitude with increasing F but the difference between 
the height of the third plateau and N{uj — 0) is changed much less. At the same time the band broadens by about 
25% (open diamonds, refer to the right vertical axis) with the uJedge given by the right hand scale in units of W/2 [see 
heavy solid line Fig. ^for a plot of N{uj) vs u) over a larger energy scale which shows the band edge]. The plateaus 
just described do not exist in an infinite band. This also holds true for the substantial temperature dependence of 
N{u} = 0) seen in the top frame of Fig. ^ as well as the thermal smearing of the phonon structure. 

The main features of renormalized (T — 0) quasiparticle density of states just described can be understood quali- 
tatively and even semi quantitatively in the context of the non selfconsistent approach. Recall (see Eq. JSJ) that the 
renormalized EDOS is 
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(29) 



We first note from Eq. ()13)) that in the clean case the imaginary part of the self energy is zero for w < wi, the first 
phonon energy in the model for a^F^oj) of Eq. 1)25)) . Hence the Lorentzian in Eq. ()29)) becomes a delta function and 
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FIG. 1: Top frame: dressed density of states N{uj) vs lj for the a'^F{uj) model of Eq. I|25^ with A = 0.71 and the bare 
bandwidth of = 2.5 eV. Impurity parameters are F = (solid) and 22.2 meV (dashed). In each group the temperatures 
are 14.5, 72.5, 145 and 435 K (from top to bottom). Bottom frame: impurity dependence of various characteristic features of 
renormalized density of states N{ijj) at low temperature T = 14.5 K. Left vertical axis: dependence of phonon plateaus height 
(starts, squares, triangles) on impurity parameter F compared with N{u!) = vs F dependence (diamonds). Right vertical axis: 

tJedge VS F. 



as a result N{uj)/N() = 1. Once u > uji but still lo < uj2 the imaginary part of the self energy becomes irai (with 
Ui = XliLLji/2). If the integral in Eq. (|29|l were not cutoff at W/2 but instead extended to infinity we would again 
get N{uj)/Nq = 1 and consequently the dressed density of states would remain unaffected by the electron-phonon 
interaction. But the finite band cutoff reduces the value of the integral in Eq. H29|l by 2ai/{W/2), i. e. by missing 
area under the Lorentzian beyond W/2. (As lo increases the integrand is no longer symmetric between positive and 
negative uj regions but we ignore this for our rough estimate so that 2ai/(W/2) is an upper limit.) As uj increases 
the Ui add until we come to the end of a^F{Lu), i. e. w = uj^ in our three delta function model. At this frequency our 
rough estimate for the reduction in N{u!) is 10% while the numerical calculations give 7.5% (see the bottom frame of 
Fig. ^ . This difference is due in part to the application of the self consistency and to our overestimate of the missing 
area under the Lorentzian of Eq. H29|l . When impurities are added, a constant uj independent term is added to the 
imaginary part of the self energy, and so in the non self consistent approximation N(lo) would now be reduced by 
2T/{W/2) at all frequencies. This expectation is in qualitative and even semi quantitative agreement with the results 
presented in the bottom frame of Fig.^ While self consistency effects are on the whole small, they are responsible for 
the fact that the lines in the bottom frame of Fig. are not quite linear in F and also not quite parallel to each other. 
At much higher energies beyond the phonon structure N{uj) drops to zero as most of the Lorenzian in the integral of 
Eq. H29|) falls outside of the range of integration. This occurs for two reasons. First, the Lorenzian becomes centered 
outside the range of integration and, second, its width becomes small. Recall that according to Eq. H13|) I]2('^ + *0+) 
is zero for u > W/2 + lod in the non selfconsistent model (pure limit). 

Finally we return to the top frame of Fig. ^ and consider more closely temperature effects. To make our main 
point it is sufficient to consider the T = 435 K curves and the value of N{lu) at o; = 0. At any finite temperature the 
imaginary part of the self energy just above the real axis is given by the expressions*^- 



5^2(0;) 
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(30) 
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FIG. 2: Top frame: the optical scattering rate Top{lu) vs uj (solid) compared with the quasiparticle scattering rate Tqp{Lu) vs 
u (dashed) for T = 22.2 meV. Temperature is T = 14.5, 72.5, 145, 425 K from top to bottom at a; = 0. The bare bandwidth 
IS W = 2.5 eV. Bottom frame: LoXopiuj) (solid) compared with negative of quasiparticle self energy — Si(w) (dashed) for the 
same parameters as in the top frame. 



which follows from Eqs. Q - ® and its w limit is 

^ ' No{0) io sinh(r!/r) 7Vo(0) ^ ^ 

The presence of this scattering rate will lead to a drop in N{uj) of Eq. H29|) by approximately 2T'^'^f^'> /{W/2), or about 
0.015 for these values of parameters, coming from the inelastic scattering which is in semiquantitative agreement with 
the numerical data. Also note the non linearity in these equations. As F increases, for example, N{llj — 0) decreases 
and thus the impurity contribution to F'^'^-'''^^ of Eq. (|31|l also decreases. As temperature is increased the inelastic 
contribution to F''^-''-^' also increases and this further reduces N{uj = 0) and consequently the effect of the impurity 
scattering. The two processes are no longer independent. 

Of primary interest in this paper is the memory function of Eqs. (|17|l - (|18|) also referred to as optical self energy. 
In the top frame of Fig. El the optical scattering rate T~p^{oj) (solid) is compared with the quasiparticle scattering 
T~p{u)) (dashed) given by — 2E2('j-'), while in the bottom frame ci;Aop(w) and —Yii{lo) are compared. The parameters 
are the same as for Fig. Q] but only results for F = 22.2 meV are presented. Many features of these curves are worth 
notice. First, the three phonon steps in the quasiparticle scattering rate (dashed), which are clearly seen in the three 
lower temperature curves, are essentially wiped out for T = 435 K (uppermost curve). For this high temperature the 
inelastic scattering due to collisions with thermally exited phonons has substantially increased the value of at 
UJ — Q above the residual scattering. It has also lead to a qualitative change in behavior at larger uj. By contrast, the 
phonon structures in the optical scattering rate at low temperature are kinks rather than steps and therefore more 
difhcult to identify. Their temperature evolution is however very similar. Complimentary to Fig. [21 in top frame of 
Fig. we compare the frequency dependence of optical (solid curves) and quasiparticle (dashed curves) scattering 
rates for several increasing values of the impurity parameter F. Shown are the results for F = 0,67, 133 meV (from 
the bottom to top) at the lowest temperature considered T = 14.5 K. In the top frame of Fig. 01 we show similar plots 
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FIG. 3: Top frame: comparison of frequency dependence of optical (solid) and quasiparticle (dashed) scattering rates at 
temperature T — 14.5 K. Impurity parameter F = 0,67,133 meV (from bottom to top at oj = 0). Middle frame: optical 
(dashed) and quasiparticle (solid) scattering rates at = vs F; dotted line shows linear dependence 27rF for reference. 
Bottom frame: Maximum value of optical (dashed) and quasiparticle (solid) scattering rates vs F; dotted line shows linear 
dependence 2n{a + F) for reference. 



for the real part of the quasiparticle self energy (dashed) compared with LoXop{uj) (solid curves). 

Note that, even for the lowest temperature considered in Fig. |21 the residual scattering in both quasiparticle and 
optical quantities are not exactly equal to their infinite band values at zero temperature, which would be 2'kT. In both 
cases it is smaller and also t~p^{llj = 0) < T^pii-o — 0). This difference between a finite and an infinite band is further 
emphasized in the middle frame of Fig. 3 where we have plotted t~j^{u) = 0) (dashed curve), t~^{lo = 0) (solid curve) 
as functions of T and compared with 27rr (dotted curve). The three curves agree in the pure limit F = but the 
deviation between these quantities increases as F increases. The order remains, with the optical scattering rate less 
than quasiparticle one, less than the infinite band value 27rF. This behavior can easily be understood from the impurity 
contribution to the imaginary part of the self energy of Eq. (|3(Jf) . As we have described before and emphasize again, 
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when r increases A'^(O) decreases so that 27rTN{0)/NQ is less than 27rr. The corresponding reduction in N{0)/No is 
roughly equal to (1 — 2r/ {W/2)) which for F — 100 meV is about 16% in good agreement with the solid curve of the 
middle frame of Fig. |31 which gives the reduction in the quasiparticle scattering rate as the impurity F is increased. 
Note that the dashed curve for the corresponding optical quantity is even lower than the quasiparticle one. This 
result comes from a full Kubo formula calculation of the conductivity and is not captured by the simplified formula 
of Eq. (|24(l for T~p^{u! = 0) which is 27rFA^(0)/A'o, the same as for quasiparticlcs. Note also that, as temperature is 
increased and inelastic processes begin to contribute to the scattering at w = 0, ''"^"'^(i^ = 0) can become smaller than 
T~p^{ui = 0) but this order is reversed as w is increased (see T — 435 K curve of Fig.|21). 

Finally we note that for higher temperatures the inelastic contribution to F^^-'^-^-' of Eq. H31|l takes the form: 



This linear in temperature law is well known and the coefficient in the square brackets would give the spectral lambda 
(A) for the infinite band case. For a finite band it is reduced as N{U)/Nq{0)) is less then one for all fi. We point out 
that as the range of a^F{^l), which is zero beyond wd, is well below (W/2) impurities and temperature will reduce 
the value of the proportionality coefficient in Eq. H32(l below its F = 0, T effective value. 

A second feature of the scattering rates shown in the upper panel of Fig.|3| which needs to be emphasized, are their 
maximum values as a function of tj. In the infinite band case both the quasiparticle and optical scattering rates would 
rise to the same asymptotic value at large uj which would be 27ra + 27rF at T = 0. This expectation is modified by the 
finite band cutoff. In the bottom frame of Fig. |21we have plotted the maximum of [T~^{u!)]max for optical (dashed) 
and quasiparticle (solid) scattering rates as functions of F and compared with 27ra + 27rF (dotted). As can be seen 
from the top frame of Fig.|21the maximum in the quasiparticle rate occurs at a frequency immediately above the third 
phonon step. For the optical rate it occurs instead at much higher values of lo. For the pure case the frequency of 
the maximum in the solid curve indeed falls beyond the bare band edge and is set not by the value of the maximum 
phonon energy, but rather by the value of the band edge itself. Also its maximum value is considerably smaller than its 
quasiparticle counterpart [by about 25%]. The deviation between the two further increases with increasing impurity 
parameter F. This difference between finite and infinite band results has important implications for the analysis of 
experimental data. Now the maximum in t~^[uj) cannot be used as a reliable estimate of the total area under the 
Eliashberg function a^F^uj), often used as a measure of the electron-phonon interaction strength. This is also the 
case for the quasiparticle rate although the differences are not as substantial. Note that the upper dashed curve in the 
top panel of Fig.|31 which gives the quasiparticle scattering rate for F = 133 meV, shows no flat region above ojd as it 
has already started to drop due to band edge effects. In fact, this is why it never reaches its infinite band maximum. 
Such band edge effects are even more substantial for the optical scattering rate which peaks only as w — > oo in the 
infinite band case. In our non selfconsistent model of Eq. (|24ll the maximum in t~^{u)) will occur at W/2 + loo- At 
this point it will have a value of approximately 27ra(l — lj\t^/{W/2)). This represents a roughly 15% reduction over 
its infinite band value in reasonable agreement with the numerical data of the lower frame of Fig. |21 

The main feature of the curves shown in Fig. |21 and |31 (top frames) that we have just described can be understood 
approximately from the non selfconsistent formulas given in the previous section. Starting with the top frame of 
Fig. the three sharp steps in T~p{uj) and the extended nearly flat region beyond are captured in Eq. (fT^ as is the 
cutoff at higher energies beyond w = W/2 + ujd with cou being maximum phonon energy (see Fig.|3| top frame), while 
the exact energy where rZ^(w) starts dropping to zero is not captured since it is due to selfconsistency that was not 
included. Similarly, Eq. H24|l allows us to understand the main differences between quasiparticle (dashed curves) and 
optical (solid curves) scattering rates. The optical scattering rate does not jump abruptly to a value of 27rai at uj ~ uji 
as T~p^{uj) does but rather grows out of zero gradually as 27rai(l — uji/lo) for frequencies in the range uji < uj < lo2- 
Additional contributions enter at uj2 and 103. On the other hand, for uj > W/2 + ujd the optical scattering rate T~p^{uj) 
does not fall off sharply but decreases towards zero as 1/uj (see Fig. O top frame). In the numerical selfconsistent 
calculations it goes faster than this and then vanishes exponentially^, but this is not captured by Eq. H24|l . 

Next, we return to the bottom panel of Fig. |21to discuss the real part of the quasiparticle self energy Si(w) [the 
negative of it is shown by dashed curves] and its optical counterpart ujXop{Lu) (solid curves). Perhaps the most striking 
feature of the real part of the self energy as a function of frequency is that it changes sign with increasing uj, as noted 
in the work of Cappelluti and Pietronerc^. In this paper we find that the corresponding memory function (or optical 
self energy) also has a "zero crossing" in a finite band. We observe that the frequency of the zero crossing ui is larger 
in the optics (about 0.5 for the parameters in this Figure) than for the quasiparticle self energy (less than 0.3). While 
the quasiparticle crossing is nearly independent of temperature in the case shown, the optical one is not, dropping 
below 0.44 at T = 425 K. 

In the top frame of Fig. ^ we show additional results for three impurity content, namely F = 0,67, 133 meV. The 
zero crossing in the memory function shifts progressively to lower frequency with increasing F, and the magnitude of 
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FIG. 4: Top frame: comparison of frequency dependence of the negative of the real part of quasiparticle self energy — Si(tj) 
(dashed curves) and corresponding memory function LjXop{uj) (solid curves) at temperature T = 14.5 K. Impurity parameter 
r = 0, 67, 133 meV (from top to bottom). Middle panel: Frequency of zero crossing a) vs F for Ei(aj) (filled boxes) and Xop{uj) 
(open boxes). Complete numerical results (solid curves) are compared with non selfconsistent estimates based on the full 
a'^F{Q,) of Eq. 12511 (dotted curves) and an appropriate Einstein oscillator spectrum (dashed curves). Bottom frame: Complete 
numerical results for optical (dashed) and quasiparticle (solid) effective mass renormalization of Eqs. 12UII and llUl vs F. Dotted 
line refers to the input A = 0.71. Dash-dotted curve shows the nonselfconsistent estimate of Eq. 1281 . equal for quasiparticle 
and optical cases. 



the maximum value of ujXop{uj) at low lu decreases correspondingly. For quasiparticles (dashed curves) the trend is 
the same. For even higher values of F than shown in Fig. ^ Aop(ti-') can become very small at small uj and even be 
negative for all (positive) frequencies. The results of our complete selfconsistent numerical calculations on the zero 
crossing frequency a) vs F dependence are summarized in the middle panel of Fig. ^ by solid curves with either open 
symbols (optics) or filled symbols (self energy) . 

The appearance of the phenomenon of "zero crossing" can be qualitatively understood with the help of the expression 
for the non selfconsistent self energy, Eq. (|12|l . and memory function, Eq. H23|l at T = 0. A crude, but reasonable 



12 



estimate can be obtained in both cases using the Einstein model for the electron-phonon spectral function: a^F(aj) = 
A S(uj — uje) with loe chosen to be a characteristic frequency of the spectrum, for example LOin- For and in the 

case wbjO) ^ W/2 (but arbitrary T and A) an explicit expression for the solution for the frequency of zero crossing 
uj can be found: 



_ I u;e{W/2) 

'"'^"V l + TM' (^^^ 

where the subscript qp means " quasiparticle" , i. e. pertinent to the self energy. It is shown in the middle panel of 
Fig.0]by dashed curve with filled boxes. For infinite band ujqp shifts to infinity and ceases to exist, as expected. Note 
that ujqp, given by this formula, become smaller as impurity parameter F increases. This behavior explains the trend 
seen in the complete numerical results (solid curve with filled boxes) . Note that the Einstein mode non selfconsistent 
estimate of the zero crossing frequency produces an underestimate of the complete numerical result. 

Similarly, for Xop{oj) we use the non selfconsistent expression Eq. H23I) with the Einstein mode a'^F{uj) and in the 
limit LUE,(jjop ^ W/2 obtain the following simple equation for Qop' 

^=2 + lni4, (34) 

qp E 

which shows that Hiop ~ V2(jjqp with logarithmic accuracy. The resulting dependence of ojop on F is shown by dashed 
curve with open boxes. Again, the general trend of the complete numerical result dependence on F (solid curve with 
open boxes) is obtained, but in this case of the optical mass renormalization the Einstein mode non selfconsistent 
estimate produces an overestimate of the exact ujop- 

The results of attempts to improve the non selfconsistent estimate for the zero crossing frequency by including the 
full elect ron-phonon spectral function a^F(S7) of Eq. I|25(l are given by dotted curves. The improvement is significant 
in the case of the self energy (dotted curve with solid boxes). This demonstrates that Cuqp depends on the shape of the 
spectrum quite strongly and is not influenced much by the self consistency. But this improvement came at a price: 
we do not have a simple formula now. In the case of the optical mass renormalization, inclusion of the full spectrum 
(dotted curve with open boxes) does not bring improvements for Uop, it even make the estimate worse as compared 
to the Einstein mode spectrum. This reminds us again of the restricted quantitative power of Eqs. H23|) and H24I) . as 
was discussed previously in Section II B. These formulae nevertheless provide a valuable qualitative guidance to the 
complete numerical results, when used properly. 

It is clear from Eqs. and (|34|1 that the zero crossing frequency contains information on the boson energy scale 
involved in the scattering process. Even though it is a slight digression from the present discussion we would like 
to point out a possible application of this finding. In their recent work on optical conductivity in high Tc cuprates 
Hwang, Timusk and GuS indeed have found a change in sign of the optical self energy (not shown in their plots). The 
frequency at which this occurs varies from compound to compound but is of the order 6000 cm~^ for their optimum 
and overdopped samples and smaller, of order 4000-5000 cm^^, in underdoppcd samples. In recent work Markiewicz— 
et al estimated that the typical band width in the oxides is of the order 1.0 to 2.0 eV for the dressed band with bare 
band structure results typically a factor two larger. This would indicate from the present work a boson exchange 
energy well above 150 meV. This is much larger than loi^ for a phonon mechanism and is consistent instead with spin 
fluctuations or marginal Fermi liquid model^''. 

Finally, turning to the position of the peaks in LuXop{u}) we note that in the infinite band case it would fall at about 
\/2lj3 for the model a^F{uj) of Eq. H25(l . Here finite band effects have shifted it down by 15% in the pure case. On 
the other hand, including F has little effect on peak's position as can be seen in Fig. 0] top frame. 

Another important characteristic of the curves in the top frame of Fig. 21 is the slope at w = 0. For the infinite 
band case it would give the mass enhancement parameter A. In the case of a finite electronic band this parameter 
can no longer be directly read off the slopes because they are changed. Denoting these by \^qj^^ and X^'ojf'^^ in the 
cases of self energy and memory function respectively (see Eq. H1U|) for example), we find that they are no longer 
equal to each other and are sensitive to impurity content, as can be seen in the top frame of Fig. ^ where we plot 
ujXop{(jj) and — Ei(a;) vs uj for three impurity content, namely F = 0, 67, 133 meV, for a low temperature T = 14.5 K. 
They also depend on temperature as shown in the bottom frame of Fig. |21 but this is not qualitatively different from 
the infinite band case. The dependence of the two X^'^^^^ on impurity parameter F for this temperature (with the 
other parameters the same as for Fig. is detailed in the bottom frame of Fig. 0] The dotted line shows the input 
A = 0.71 for reference. Both optical (dashed) and quasiparticle (solid) mass renormalization decrease substantially 
with increasing F and the former quantity is always smaller. Note that in the non selfconsistent case Eq. H14|l applies 

to both X^qJ^'^ and X^Jjf'^K This dependence is also shown by long dash - dotted curve for comparison. While in the 
pure case it agrees well with the exact result for the quasiparticles, it deviates substantially from the exact optical 
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FIG. 5; Top frame: real part of optical conductivity ai{uj) vs lu for F — 44, 67, 89, 111, 133, 222 meV (from top to bottom at 
u — 0). Middle frame: dependence of the optical sum S of Eq. 1161 on 27rF (solid) and on To,^{lo = 0) (dashed). Bottom frame: 
dependencie of dc resistivity p — l/ai{Lo — 0) on 27rF (solid) and on Top{Lo = 0) (dashed). Dotted curve is the infinite band 
result. Temperature is T = 14.5 K for all frames. 



curve. As T increases the deviations increase although the trend is properly given. This difference described goes 
beyond the approximations that we used to obtain Eq. (|28() from the exact expression of Eq. based on the Kubo 
formula for the optical conductivity. These approximations are clearly not very accurate but because the resulting 
formulas are quite simple and analytic they are nevertheless useful. 

In the top frame of Fig. |31 we show results for the real part of the conductivity ai{uj) vs lo at temperature T ~ 
14.5 K. Optical conductivity is measured in units ire^ /\2Dm{W/2)]. Six values of impurity parameters are used, 
namely F = 44.4, 66.6, 88.8, 111.0, 133.1, 221.9 meV. There are several features of these curves which are different from 
corresponding infinite band results. Perhaps the most obvious is that the total optical spectral weight S, defined in 
Eq. H16() . is no longer independent of temperature and impurity content. [Temperature dependence is discussed in our 
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previous paper—.] The dependence of S* on F is shown by the solid curve in the middle frame of Fig. |S| The optical 
spectral weight is measured in units ire^ l[2Dm], such that in the infinite band case S — 1/2. Because residual and 
inelastic scattering is no longer strictly additive in finite bands, the impurity parameter F cannot be directly obtained 
from optical conductivity experiment. Therefore, we also plot S vs t~j^{uj = 0) (dashed curve), a parameter that is 
measurable. While there is a small difference between the two plots, both show that the total optical spectral weight 
is reduced as F is increased. Another important property of the conductivity cri(w) from Fig. top frame, which 
needs to be commented upon, is its dc value, or its inverse l/cri(w = 0), the resistivity. This quantity is plotted in 
the bottom frame of Fig. O where it is seen to increase with F (solid curve) . This is also the case when plotted with 
respect to t~j^ (oj — 0) (dashed curve) , and the dependence is not quite linear as would be expected in an infinite band 
(dotted hue). 

The behavior of solid curve can be understood from the formula for the resistivity: p{T) = [tt/{2S)]t^p^{uj = 0). 
While, as we have seen in a previous section, the approximate analytic formulas that can be obtained for the optical 
quantities, are not as accurate as for the self energy nevertheless they can be quite helpful in providing insight into 
the complete numerical results. Sharapov and Carbotta^l give the following expressions for inelastic and impurity 
contributions, respectively: 

V!pho„(^-0) = dna'Fm J^^de^^[nin) + fin-e)] (-^) , (35) 

-op,i„.p(- -> 0) = 2nTj_^ de^^ ' ^^^^ 

Besides the explicit thermal factor appearing in these equations, the temperature also enters through the renormalized 
DOS factor N{e) which also depends on the phonon a^F(il) and on F in contrast to the infinite band case. Independent 
of the details, because the N{e) factor is everywhere smaller than its infinite band value of one, the resistivity is always 
reduced below its infinite band value. This reduction increase with increasing value of F as we see in lower frame of 
Fig.Oand is also increased with increasing temperature. At low temperatures the appropriate measure of the decrease 
in the impurity term of Eq. (|36|l is the value of N{0) while for the inelastic term it depends on N{^1) with fl within 
the phonon range. 



IV. VERY NARROW BANDS 



The case considered so far corresponds to a rather broad band as compared with the phonon energy. Nevertheless 
we found important qualitative changes from the infinite band case. Band structure calculations^i for K^Cqo, as 
an example, give a half band width of about 250 meV which is now comparable to the energy of the maximum 
phonon energy of 190 meV in our model electron-phonon spectral density. For such cases Kostur and Mitrovio^ 
and later Pietronero, Strassler and Grimaldi^S have considered the effect of vertex corrections and a generalization 
of the Eliashberg equations which go beyond the Midgal theorem. The specific case of the Pauli susceptibility was 
considered by Cappelluti, Grimaldi and Pietronero42.. In more recent work, Cappelluti and Pietronero^ recognized 
that it was the effect of a finite band that primarily accounted for some of the qualitative differences found in their 
previous work and proceeded to include only these as a first step in understanding self energy renormalization. Here 
we follow their lead, but consider instead optical properties. 

In this section we wish to accomplish three goals. First, we want to understand differences that arise when very 
narrow bands are involved as compared with relatively wide ones. Second, we want to compare a case with a larger 
value of A and finally we replace the three delta function model for a'^F{uj) with a more realistic extended spectrum. 
We consider a model with three truncated Lorentzians: 



^ ' ZTT 

i=l 



e(?7,-|a;-a;,|), (37) 



where Ui, the centers of the peaks, are the same as in Eq. (|25|l . For each peak the parameter Si controls the half 
width, while the full spread is equal to 2ryi. The rescaling factor i?(A) is inserted to guarantee a chosen value of A. 
Truncated Lorentzians are often used in the literature to introduce a smearing of the simple Einstein mode spectrum. 
For our numerical work we picked Si = 0.2u;i and r]i = 0.611;^. In this case the peaks are wide and overlapping. The 
spectrum of Eq. (|37() has the characteristic logarithmic^ frequency of tj;„ = 96 meV. 

While, to set the parameters used here, we consider what might be reasonable for K^Cqq, i. e. we chose W/2 — 
250 meV and an effective A of about one, we do not imply that our calculations can be applied directly to this specific 
system. Other complications such as the effect of Coulomb interactions-i may need to be included as well. For 
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FIG. 6: Frequency dependence of the renormalized density of states N{uj)/No for a very narrow band with W/2 = 250 meV 
(thin curves) with the parameters described in detail in the text. The heavy curve refers to N{uj)/No with W/2 = 1.25 eV for 
comparison. 

example the small Drude peak seen in the experiments of Degiorgi et aP^^ which has a width of about 10 — 20 meV 
and a weight representing only about 10 — 20% of the total spectral weight is not understood in our work. On the 
other hand qualitative features of the memory function vs frequency dependence, such as the observed zero in its real 
part, are captured. 

In Fig. El we present a series of results for the dressed quasiparticle density of states N{uj)/No as a function of 
energy cu. The heavy solid curve shows previous results for a rather wide band W/2 = 1.25 eV on a broader scale for 
the three delta function model of Eq. H25|) with A = 0.71 and T = 22 meV at low temperature T = 1.25 K. It is to be 
compared with the other curves all of which are for a band that is five times narrower, namely W/2 = 250 meV. While 
for the wider band significant phonon structures are limited to an energy region well below the bare band cutoff at 
uj/{W/2) = 1, for the narrower band they dominate the shape of N{u})/No even beyond lu/{W/2) = 2. No particular 
signature associated with the bare band edge remains. This is distinct from the heavy continuous curve which shows 
a smooth drop off at a new easily distinguishable renormalized band edge energy increased somewhat over its bare 
value and smeared by the interactions. 

All thin solid curves in Fig. |B|are for the three delta function model of a'^F{uj) with A = 0.71 at low temperature 
(T = 1.25 K). The impurity parameter is F = 4.4,13.3 and 39.9 meV (from top to bottom). The dashed curve 
correspond to the extended spectrum of Eq. H37|) with A — 0.71 and 7 = 5 meV at low temperature (T = 2.15 K). The 
two sharp step like drops at ijj/{W/2) < 0.40 present in the delta function case are now almost completely smeared 
out. The sharp spike like minimum at the energy of the maximum phonon energy u!/{W/2) = 0.76 [to = 190 meV] 
and the near vertical drop at twice this energy, seen in the solid curves are gone in the dashed curve as are the distinct 
multiphonon structure at higher energies. The plateau like region at 0.8 < uj/{W/2) < 1.3 [for oj between 190 and 
390 meV] in the solid curves becomes a shoulder in the dashed curve. 

The modification of the renormalized density of states depends on the mass renormalization parameter A. To 
demonstrate this we present in Fig. Elthe result for the same extended spectrum of Eq. (|37|l ) and the same impurity 
parameter F = 4.4 meV but with A = 2 (the dotted curve). Now the low temperature step at uj/(W/2) = 0.16 
becomes visible and a deep and wide minimum develops in the frequency region of the most strongly coupled part of 
the electron-phonon spectral density centered at lu/{W/2) = 0.76. Additionally, the multiphonon processes become 
stronger, which is manifested by appearance of the maximum at ll!/{W/2) « 1.4 in the dotted curve in place of the 
shoulder in the dashed curve. Finally, the renormalized band edge has been shifted to much higher frequency and is 
not visible in Fig. |S| 

Returning to the light continuous curves we note that, compared to the wide band case of the previous section, the 
phonon steps are now significantly reduced in magnitude as the impurity scattering is increased. For example, compare 
the bottom curve for F = 39.9 meV to the top one for F = 4.4 meV. These reductions go beyond non selfconsistent 
approximation and demonstrate that the selfconsistency becomes more important as the bare band width is reduced. 
The near additivity of electron-phonon and impurity contributions is lost. While for the curve corresponding to the 
purest case, the non self consistent approximation predicts well the size of the first step, it is not as good for the 
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FIG. 7: Top panel: Comparison of optical (solid curves) and quasiparticle (dotted curves) scattering rates. Bottom panel; 
Comparison of minus the real part of the self energy (dotted) with the corresponding optical quantity uXopiui) (solid). The mass 
renormalization is A = 2.0 (grey curves) and A — 0.71 (black curves). In all cases there are two sets of curves corresponding to 
impurity parameter F = 5 and 25 meV. The extended a^F{Q) of Eq. 13711 was used. The half band width W/2 — 250 meV, 
temperature T = 2.15 K. 



second and the third step is quite off. This is expected as in this case we are already not so far from the bare band 
edge and the selfconsistency becomes essential. 

Note that at w = for the bottom light soUd Une with F = 39.9 meV, N{0)/Nq sa 0.77 in the numerical works. 
The non selfconsistent estimate (1 — 2r/{W/2)) « 0.68 which is considerable smaller. However self consistency has 
actually reduced the impurity scattering rate below its infinite band value of 2ttT because it is equal to 27rF7V(0) and 
A^(0)/A'o(0) is smaller than one. Accounting for its 23% reduction eliminates much of the difference described above. 
A similar semiquantitative argument can be made to understand the reduction in phonon step size as F increases 
when the simple non selfconsistent estimate begin to fail. Finally, we note the crossing of the light solid curves around 
Lo/{W/2) sa 0.8 and the increase in density of states in the tails beyond uj/{W/2) — 1.5 as the impurity scattering is 
increased. 

In Fig.[7|we turn to the memory function (solid lines) which is compared with the quasiparticle self energy (dotted 
lines) for several cases. We present results for A = 0.71 (black curves) and for A = 2.0 (grey curves) and for two 
different impurity content, F = 5 and 25 meV calculated with the model of the extended a^F{fl) of Eq. H37() . As 
expected both quasiparticle and optical scattering rates rise to a higher maximum value when A is larger, but the 
increase is not linear. This is followed by a drop towards zero as to gets large instead of saturating at a common value 
of 27ra + 27rF as discussed previously. The intercept of t~p^{lo) and t~p^{lj) at cj = is related to the elastic impurity 
scattering as the temperature for the figure is small T = 2.15 K. For T~p^ it does not depend on the value of A but 
for T~p^ it does. This can be seen most clearly in the top set of curves for which F = 25 meV but is noticeable in the 
lower set with F = 5 meV. When A is larger T~p^ at cj = is smaller. This effect results from the Kubo formula and 
is not captured in any of our approximate analytic formulas. As previously noted T~p{uj = 0) < T~p(ijj = 0) < 27rF. 

For A = 2 the maximum value of T~p{bj) (see grey dotted curves in the top panel of Fig. 0) occurs at a frequency 
slightly below u)/{W/2) — 1 and is almost the same for the F = 5 meV and F = 25 meV. Although we have increased 
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the impurity scattering we have not gained in maximum quasiparticle scattering. In an infinite band it would have 
risen from 3.22 to 3.72 for the parameters used. This effect is due to the self consistency. Both elastic and inelastic 
scattering involve not just a^F{U) and F respectively but also the self consistent value of the dressed quasiparticle 
density of states N{u!)/Nq which becomes reduced as F is increased. This in turn reduces both elastic and inelastic 
scattering rates leading, in case considered, to a saturation of the maximum in T^p^{uj). For F = 5 meV the rate is 
about 34% lower than its infinite band value and for F = 25 meV it is about 43% lower. A similar situation holds 
for the case when A = 0.71 (black dotted curves) although in that case the maximum quasiparticle scattering does 
increase slightly with increasing F, for F = 5 meV it is 26% below its infinite band value and for F = 25 meV it is 40% 
below. Turning next to the optical scattering rate (solid curves in Fig. Q top frame), we note first that they peaks at 
a higher frequency than does the corresponding quasiparticle rates. Also they are considerably smaller in magnitude, 
approximately 0.95 and 0.8 for F = 25 and 5 mcV respectively for the case A = 2. 

In the bottom frame of Fig.[7|we compare our results for the minus real part of the electronic self energy — Si((jj) 
(dotted curves) with the corresponding optical quantity ujXop{uj) (solid curves) of Eq. H18|l. which is the minus of the 
imaginary part of the memory function M(w) defined in Eq. I|15(l . It is quite clear that the optical masses (slopes at 
u! = for solid curves) are in all cases considerably smaller than the quasiparticle masses (slopes at a; = for dotted 
curves). As an example, for A = 2 the quasiparticle mass at F = 5 meV is 1.1 as compared to 0.70 for the optical 
mass, while at F = 25 meV we have 0.75 and 0.35 respectively. Using the non self consistent formula Eq. (|28|l gives 
for both masses A = 2 (A = 0.71) 1.1 (0.39) and 0.94 (0.21) for F = 5 meV and F = 25 meV respectively. For the 
pure case the agreement for the quasiparticle mass is good but this is no longer the case for F = 25 meV. Also for 
the purer case considered on Fig. the zero crossing of quasiparticle and optical mass renormalization function can 
be understood qualitatively with Eqs. (|12|l and H23|) but these simple estimates begin to fail for higher values of F. 

Finally, we comment on the reflectivity data of Degiorgi et a^^^ on K^Cqq. They did not analyze their data 
to extract optical scattering rate t~p^{lu) amd mass renormalization uiXop{uj). Nevertheless we infer from the data 
presented three qualitative features. The value of rj^^(w) at w = which gives a measure of the residual scattering is of 
order 160 meV. This large value is incompatible with the observed small Drude peak in (Ti(a;) of width 20 meV which 
contains about 12 % of the total optical spectral weight. Second, at a; ~ 500 meV the scattering rate has increased 
to approximately 500 mcV which implies an inelastic contribution of 360 meV. Such a rise is much larger than can be 
achieved in the model of Fig. 6 and would indicate that the bare band width is somewhat larger than present band 
structure calculations predict and that the spectral A defined by the input electron-phonon spectral function a'^F{Lu) 
is even larger than 2. Thirdly, luXop{ll}) changes sign at approximately 220 meV after which it plunges towards large 
negative values. This feature can be taken as the hallmark of finite band effects as it does not occur in infinite band 
theories. While such a zero crossing occurs naturally in our calculations and is generic, it is not clear that a set of 
microscopic parameters chosen to reproduce the features of the scattering rate T^p^{uj) would also accurately produce 
the position of the zero in the optical mass. We did not attempt such a combined fit as it would require a value of A 
which appears to be rather large. 

V. CONCLUSIONS 

In an electronic system with a constant bare electronic density of states A*o the application of a finite band cut off 
profoundly modifies the electron-phonon renormalization effects. In the infinite band case the dressed quasiparticle 
density of states N{u;) remains equal to A^o and is independent of impurity scattering and temperature. A self 
consistent solution for the self energy in a finite band shows that N{lu) acquires low energy structure on the scale of 
the phonon energies. The band edge becomes smeared and the band extends beyond the original bare cut off. This 
extension of the band to higher energies is accompanied by a compensating reduction of spectral weight below the 
bare cut off. Equally importantly N{uj) is affected by impurity scattering and by temperature. N{u! = 0) is reduced 
in both cases. On the other hand, while temperature rapidly smears out the phonon structure, impurities mainly 
reduce its amplitude. 

The emphasis of previous works was on the effects of temperature and impurity scattering on the electron self 
energy S(w), which is the quantity that determines quasiparticle properties. Here we have extended these works to 
optical properties and considered characteristic features of the memory function. For the infinite band case elastic 
impurity scattering just adds a constant amount (27rF) to the quasiparticle inelastic scattering rate, but in our case, 
because of the application of self consistency, they no longer add. Even at = we find that T~p < t^^ < 27rF. The 
well known result that, at high energies both optical and quasiparticle scattering rates due to phonons become equal 
and saturate at a value 27ra (with a the area under the electron-phonon spectral density) no longer holds. While the 
maximum in can come close in value to 27ra + 27rF, the corresponding optical quantity is much smaller in 
magnitude. Its maximum value increases with increasing a but this increase is sublinear. A similar situation holds 
when F is increased. At yet higher energies both scattering rates go to zero because of the finite band. 
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The known result that the real part of the quasiparticle self energy is unaffected by the impurity scattering no longer 
holds and this has an impact as well on the optical mass renormalization. Quasiparticle and optical effective mass 
renormalization at w = now differ from each other, neither is equal to A and both depend on impurity scattering. 
The real part of the self energy changes sign with increasing energy as does a;Aop(ct'). The energy at which the zero 
crossing occurs is set by the phonon energy scale and can depend both on temperature and impurity content. It is 
larger for optics than it is for the self energy in many cases but not always. The optical spectral weight, i. e. the 
area under the absorptive part of the optical conductivity varies with temperature and with impurity scattering. The 
elastic and inelastic contributions to dc resistivity are no longer additive. 

All these effects were found to be significant in magnitude even when a rather modest value of mass enhancement 
parameter A = 0.71 is used with a band width of 2.5 eV. While a three delta function a^F{uj) was used to emphasize 
boson structure with maximum phonon energy of 190 meV, a broader spectrum was also considered. This soften 
phonon structures but did not eliminate them. Of course, for simple metals such as Pb, the band width is much wider 
than considered above and the maximum phonon energy is also an order of magnitude smaller, so that in this case 
the infinite band approximation is appropriate and the finite band corrections found in this paper would be negligible. 
We have found that increasing the value of A to 2 increases boson structure but the increase is not linear in the value 
of A. Also decreasing the value of W to 500 meV, a value suggested by band structure calculations in the alkali doped 
Cgo, leads to a new regime in which important modifications due to the electron-phonon interaction dominate at all 
energies and no easily identifiable trace of the underlying bare electronic band cutoff remains. 

While all these conclusions are based on numerical solution of the selfconsistent equations for the self energy and 
the Kubo formula for the conductivity, we have also derived more transparent analytical formulas evaluated in a 
non selfconsistent approximation. These simple formulas are not always accurate but give considerable insight into 
complete numerical results obtained and prove valuable in the analysis of optical data in finite band metals. 
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